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CHEMISTRY.  KINETICS  AND  THERMODYNAMICS 
James  H.  Bigelow* 

The  Rand  Corporation,  Santa  Monica,  California 


I .  INTRODUCTION 

In  this  paper,  the  first  of  a  series  of  three  ,  we  will 
discuss  the  equivalence  of  three  ways  of  viewing  the  classi¬ 
cal  chemical  equilibrium  problem*  From  the  point  of  view 
of  classical  chemistry,  this  problem  can  be  stated  as  follows: 
Given  the  set  of  reactions  that  may  occur  among  all  the 
species  of  a  chemical  system,  their  equilibrium  constants, 
and  the  initial  composition  of  the  system  (i.e.  the  amounts 
of  each  species  present  in  the  system  initially),  find  the 
equilibrium  composition  of  the  system. 

From  the  point  of  view  of  kinetics,  the  problem  is  to 
determine  the  evolution  of  the  various  parameters  describing 
the  system  from  their  initial  values.  When  the  parameters 
cease  to  change,  the  system  must  be  in  equilibrium. 


1  y?  1,r  ' 1  r  ■ 

Any  views  expressed  in  this  paper  are  those  of  the 
author.  They  should  not  be  interpreted  as  reflecting  the 
views  of  The  Rand  Corporation  or  the  official  opinion  or 
policy  of  any  of  its  governmental  or  private  research  sponsors. 
Papers  are  reproduced  by  The  Rand  Corporation  as  a  courtesy 
to  members  of  its  staff. 

‘The  other  two  papers  are,  in  order,  "Degeneracy  in 
Ideal  Chemical  Equilibrium  Problems,"  and  "Computing  Equilib¬ 
rium  Compositions  of  Ideal  Chemical  Systems." 
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Given  the  same  system,  thermodynamics  seeks  that 
composition  at  which  the  usable  (or  free)  energy  remaining 
in  the  system  is  at  a  minimum,  It  can  be  shown  that  the 
system  will  then  be  at  equilibrium. 

Much  of  the  materiel  contained  herein  is  old.  It  is 
presented  here  in  a  consistent  notation,  and  in  a  way 
intended  to  emphasize  the  common  aspects  of  the  three  points 
of  view. 
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II.  CLASSICAL  CU KM IS TRY 


Classically,  the  composition  of  a  single-  or  multi¬ 
phase  chemical  system  in  equilibrium  has  been  computationally 
determined  by  solving  certain  systems  of  simultaneous  equa¬ 
tions.  These  equations  included  linear  equations  (mass 
balance  laws)  as  well  as  nonlinear  equations  (mass  action 
laws).  It  is  possible,  however,  to  formulate  the  problem 
as  one  of  minimizing  a  particular  nonlinear  function  (the 
free  energy)  subject  to  linear  constraints.  This  was  first 
done  by  J.  Willard  Gibbs  [l],  in  his  famous  paper  "On  the 
Equilibrium  of  Heterogeneous  Substances."  It  was  later 
shown  in  a  different  way,  and  the  equivalence  proved  with 
greater  mathematical  rigor  by  Shapiro  and  Shapley  [2], 

This  section  restates  the  first  part  of  Shapiro  and  Shapley. 

1.  The  Chemical  System 

Vie  will  consider  chemical  systems  composed  of  a  finite 
number  (not  excluding  one)  of  homogeneous  phases.  A  ho¬ 
mogeneous  phase  is  a  mixture  of  chemical  species  which 
is  homogeneous  in  chemical  composition,  pressure,  and 
temperature.  The  question  of  under  what  circumstances 
a  phase  should  be  regarded  as  homogeneous  can  only  be 
answered  as  part  of  the  process  by  which  we  formulate  a 
model  of  an  actual  chemical  system. 

For  example,  a  chemical  system  consisting  of  a  vapor 
over  a  liquid  solution  might  be  regarded  as  having  two 
phases,  a  gas  phase  and  a  liquid  phase.  A  chemical  system 


consisting  of  two  solutions  separated  by  a  semipermeable 
membrane  might  be  regarded  as  having  two  liquid  phases. 

All  this  model  requires  of  a  phase  is  homogeneity.  A 
phase  need  not,  for  example,  occupy  contiguous  portions  of 
space.  Thus  in  [3],  the  interiors  of  all  the  red  blood  cells 
of  the  body  are  usefully  regarded  as  forming  a  single  phase. 

The  entities  in  a  phase  are  called  species,  by  which 
we  mean  molecular  species.  Examples  of  species  are  H20, 

NaCl,  Na**\  When  a  species  is  defined,  its  molecular  or 
ionic  structure  is  implied,  as  well  at  the  number  of  atoms 
of  each  element  composing  a  molecule  of  the  species.  Thus 
the  molecular  formulas  of  two  distinct  species  might  be 
identical.  (An  example  of  this  would  be  right-handed  and 
left-handed  amino  acids  . )  Also,  we  will  find  it  mathemat¬ 
ically  convenient  to  regard  a  species  which  can  occur  in 
two  different  phases  as  two  different  species.  Thus  if 
H£0  can  occur  in  both  a  liquid  phase  and  a  gaseous  phase, 
we  would  refer  to  two  different  species  H20,  perhaps  calling 
them  H20  liquid,  and  H20-vapor.  When  a  molecule  moves  from 
one  phase  to  another  (by  evaporation,  condensation,  migra¬ 
tion  across  a  semipermeable  membrane,  or  black  magic)  we 
will  regard  a  type  of  chemical  reaction  as  having  occurred. 

By  saying  that  a  species  can  occur  in  a  given  phase, 
we  do  not  mean  that  at  equilibrium  any  positive  amount  of 
that  species  will  be  found  in  the  phase.  We  only  mean 
that  wo  have  made  provision  for  the  possibility.  The 
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question  of:  which  species  should  be  provided  for  in  what 
phases  can  only  be  answered  by  judgment  and  experience, 
and  by  some  knowledge  of  the  system  under  study. 

It  should  be  pointed  out  that  if  a  phase  is  a  chemical 
solution,  then  the  solvent  itself  is  one  of  the  species 
occurring  in  that  phase. 

We  number  the  species  in  the  chemical  system  1,  2, 

. . . ,  n.  Let  be  the  number  of  moles  of  species  j  in  the 
system.  (By  definition,  a  mole  is  Avogadro's  number  of 
molecules,  approximately  6  x  10  .)  We  denote  the  phase 

containing  species  j  by  <j>,  and  we  may  indicate  that 
species  j  and  k  are  in  the  same  phase  by  writing: 

j  e  <k>,  or  k  «  <j>  or  <j>  »  <k>. 

Each  phase  has  associated  with  it  a  sum. 

(II.  l.D  x<,  o  S  x,. 

K  j  g  <k>  J 

Each  species  has  associated  with  it  a  mole  fraction. 
(II.  1.2)  x  » 

It  is  easy  to  see  that  x4  is  the  concentration  (on 

J  " 

the  mole  fraction  scale)  of  species  j  when  the  system  has 
composition  x  «  (x^,  x2  ,  .  . .  ,  x  ) . 


2.  The  Mass  Balance  Laws 


Following  [4],  we  will  express  the  mass  balance  laws 
as  linear  conditions  on  x.  Let  B^,  . ..,  B^,  .  ..*  Bm  be 
a  set  of  fundamental  building  blocksj  such  that  each  species 
j  is  a  unique  combination  of  these  building  blocks. 

For  example,  B^,  . ..,  B^,  . ..,  Bm  might  be  taken  to 
be  the  atomic  elements.  However,  it  is  often  better  to 
choose  more  complex  structures  as  building  blocks.  (See 
[4]  and  [5].)  Let  each  molecule  of  species  j  contain  a^ 
units  of  B^.  Thus,  the  vector  (a^,  *2 y  *  *  *  *  araj) 
essentially  the  molecular  formula  for  species  j.  Note 
there  is  nothing  to  prevent  two  different  species  from 
having  identical  molecular  formulas.  Let  b^  be  the  total 
number  of  units  of  B^  in  the  system.  We  then  have  the 
conditions : 


(XX.  2.1) 


n 


A  auxj  ■  bt- 


The  description  of  those,  aspects  of  the  system  con¬ 
cerned  with  such  questions  as  which  species  are  permeable 
to  interphaso  boundaries  can  be  implemented  by  an  appro¬ 
priate  extension  of  the  set  of  fundamental  building  blocks. 
This  process  has  been  described  elsewhere  [4] .  In  addi¬ 
tion,  it  may  be  desirable  to  impose  a  condition  of  elec¬ 
trical  neutrality  on  one  or  more  phases.  This  too  may  be 
expressed  in  the  form  of  (II. 2.1)  ([2]  p.  357).  Other 
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constraints,  such  as  thGse  on  the  volume  of  a  particular 
phase,  can  be  similarly  handled. 

Finally,  from  the  definition  of  x^,  it  is  clear  that 


i<re  must  have 


3.  Reaction  Vectors  and  the , Stoichiometric  Conditions 
In  equilibrium  chemistry,  a  reaction  is  described 
'by  writing  two  formal  sums  separated  by  a  double  arrow* 

If  we  denote  by  Nj  the  name  or  chemical  symbol  of  species 
j ,  then  a  reaction  would  be  described  by  thfi  "equation" 

,,(11.3.1)  E  r.  K.  E  p.  N.  . 

"  J  J  J  J 

Each  of  the  terms  in  (II. 3.1)  consists  formally  of  a  posi¬ 
tive  rev  number,  r ^  or  p  ^ ,  called  a  stoichiometric  coef— 
f ip lent,  multiplied  by  the  symbol  for  some  chemical  species 
of  the  system.  The  species  occurring  in  the  left-hand 

formal  sum  with  a  positive  coeefficient  r.  are  called 

3 

reactants;  the  species  occurring  in  the  right-hand  formal 
sum  with  a  positive  coefficient  p^  are  called  products . 

The  two  formal  sums  separated-  by  a  doub3:c  arrow  form  the 
stoichiometric  equation  of  the  reaction.  Note  that  math¬ 
ematically  speaking,  the  stoichiometric  equation  is  not 
an  equation  at  all,  unless  the  l\L  are  interpreted  as  vectors. 

We  can  measure  the  extent  to  which  the  reaction  (II. 3.1) 
takes  place  by  an  extent  of  reaction  parameter  §.  "Equation'* 
(II. 3.1)  is  interpreted  to  mean  that  if  the  extent  of  renew, 
tion  g  changes  to  g  +  d§,  then  an  amount  r-^dg  of  species 
Np  r2^  sPGciGS  •••j  and  rnd§  species  N  have 
reacted  to  form  an  amount  p-^dg  of  species  Np  p2dg  of  species 
N«,  ...,  and  p  dg  of  species  N  . 

^  li  n 
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Given  any  such  stoichiometric  equation,  we  shall  de¬ 
fine  a  vector  0,  with  components  0.,  (1  £  j  £  n)  by  viewing 

J 

the  expression  (II. 3.1)  as  though  it  were  a  mathematical  vec¬ 
tor  equation*  and  treating  the  double  arrow  (♦*)  as  though  it 
were  an  equals  sign  (»)'..  Bring  all  the  terms  to  the  left 
(changing  the  signs  in  the  process),  and  collect  the  coefficients 

of  each  N..  Then  if  we  let: 

3 

8i '  rj  -  pj' 

we  will  have  the  "equation": 

(11. 3. 2)  S  0,  N.  *  0. 

j=l  J  J 

We  call  any  vector  0  so  obtainable  from  a  reaction  a 
reaction  vector.  Observe  that  different  reactions  may 
have  the  same  reaction  vector. 

As  a  condition  the t  0  be  a  reaction  vector,  we  have 

that: 

n 

(11. 3. 3)  E  a,,  e,  «  0  (1  <  i  <  m) 

3=1  1J  J  “  ~ 

where  the  a...  arc  those  of  equations  (II. 2.1).  The  condi- 

J 

tion  can  be  made  sufficient  as  well  as  necessary;  in  actual 

practice,  the  model  maker  judiciously  designs  (II. 2.1) — 

that  is,  determines  the  a.. — so  that  (11,3.3)  will  be  a 

•*•3 

sufficient  as  well  as  a  necessary  condition.  (For  a  simple 
example,  see  [2],  pp.  338-359.)  We  shall  regard  any  vector 
0  that  satisfies  (11.3.3)  as  a  reaction  vectoe, 

Q 


Mhite  (II. 2.1)  and  (II. 3. 3)  are  closely  related,  they 
are  not  identical.  Unfortunately,  the  chemical  literature 
refers  to  both  sets  of  conditions  as  mass  balance  laws. 

We  shall  reserve  the  term  'faass  balance  laws"  for  (II, 2.1) 
and  call  (II. 3. 3)  the  stoichiometric,  conditions. 
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4,.  The  Mass  Action  7, awn 

For  the  time  being,  we  shall  confine  our  attention 
to  systems  which  exhibit  the  simplest  form  of  the  mass- 
action  -laws.  Such  systems  are  called  "ideal"  in  [ 6 ] .  A 
further  discussion  of  ideality  can  be  found  in  [7-3..  We. 
will  discuss  nonideal  systems  in  a  later  section. 

The  mass  action  laws  can  be  stated  as  follows:  For 
any  reaction  having  a  reaction  vector  0,  there  is  an  equi¬ 
librium  constant  lc(6),  such  that,  for  any  composition  vec¬ 
tor  x  representing  an  equilibrium  state  of  the  system,  we 
have  ■ 


(II. 4 .1) 


n  0, 

n  i,  -1  »  k(0). 


Note  that  k(o)  does  not  depend  on  x.  (II. 4.1)  is  a  condi¬ 
tion  that  x  must  satisfy  if  it  is  to  be  an  equilibrium  com¬ 


position. 


Wo  shall  find  it  useful  to  perform  certain  manipula¬ 
tions  on  (11.4,1).  Leu  us  first  take  the  logarithm  of 
each  side  to  get : 

n 

(II. 4. 2)  X  e.  log  x.  «  log  k(e), 

j-1  J  J 


Equation  (II. 4. 2)  is  thus  satisfied  for  all  reaction 
vectors  q.,  that  is,  for  all  vectors  0  which  satisfy  the 
stoichiometric  conditions  (II. 3. 3).  It  should  be  noted 
that  the  form  of  (II. 3. 3)  implies  that  if  a  and  m  are  two 
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reaction  vectors,  then  (ae  +  Pep)  is  also  a  reaction  vector, 

for  any  real  numbers  a  and  P.  By  the  same  token,  the  form 

of.  (II. 4.2)  shows  that  for  the  same  0  and  cp, 

log  k(ae-  +  pep)  »  a  log.  k(e)  +  p  log  k(cp) . 

Thus  log  k(e)  is  a  linear  function  of  0,  so  that  there 
mur.i  exist  constants  c^  such  that: 

n 

(11. 4. 3)  log  k(G)  «  -  E  c,0.. 

j«l  J  J 

Substituting  the  value  of  log  lc(o)  from  (II. 4. 3)  into 
(II. 4. 2)  allows  us  to  restate  the  mass  action  laws  as  fol¬ 
lows:  There  are  constants  c.  (1  <  j  <  n)  such  that  if  0 

J 

is  any  reaction  vector,  and  if  x  is  any  equilibrium  com¬ 
position,  then 

n 

(11. 4. 4)  E  0.(c.  +  log  x.)  =  0. 

j«l  J  J  3 

This  form  of  the  mass  action  lav?,  which  is  quite  equiva¬ 
lent  to  the  standard  forms,  is  more  convenient  for  our 
purpose. 

Another  form  is  often  useful  for  computational  pur¬ 
poses.  Consider  (c  +  log  A)  as  an  n-dimensional  vector, 
whose  j —  component  is  c^  +  log  k  .  If  we  also  think  of 
the  rows  ai  of  the  matrix  A  as  n-vectors,  then  (II. 4. 4)  may 
be  restated  as:  If  x  is  an  equilibrium  composition,  then 
(c  +  log  x)  is  orthogonal  to  every  vector  n  which  is  in 
turn  orthogonal  to  all  the  vectors  a.^.  But  it  follows 
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frbin  elementary  linear  algebra  that  thin  condition  is  sat¬ 
isfied  if  and  only  if  (c  +  log  x)  is  representable  as  a 
linear  combination  of  the  rows  a^  of  A.  That  is,  con¬ 
dition  (II. 4. 4)  is  satisfied  if  and  only  if  there  exists 
a  vector  tt  ■  (rr^,  tt •••*  trm)  such  that,;  in  matrix  nota¬ 
tion, 

(II. 4.5)  c  +  log  x  »  A^tt. 
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5.,,  The  Gibbs  Free  Energy  Function 

We  will,  have,  frequent  occasion  to  refer  to  the  free 
energy  of  a  system.  In  Section  IV  we  will  discuss  this 
more  fully;  for  now,  we  note  that  since  we  are  dealing 
with  ideal  systems  with  each:  phase  at  constant  temperature 
and  pressure,  the  free  energy; may  be  written  as  a  function 
of  composition  as  follows  [6],  [.7]' 

n 

(II. 5.1)  F(x)  *  £  k,(c.  *  log  £.). 

j=*l  J  J  v 

The  numbers  c.,  are  called  free;  energy  parameters,  but  still 
they  satisfy  (11.4.3)* 

The  least-action  principle  ([6l,  p,  29)  for  the  free 
energy  states  that  the  system  is  in  equilibrium  if  and  only 
if  its  free  energy  is  a  minimum,  subject  to  the  constraints 
(II. 2. 1-2) ;  that  Is,  a  composition  vector  x°  represents  an 
equilibrium,  state  if  and  only  if  F(x°)  <  F(x)  holds  for  all 
states  x  which  satisfy  the  conditions  (II. 2.1)  and  (II. ,2. 2) 
and  which  are,  sufficiently  close  to  x°. 

Note  that  the  least-action  principle,  as  stated  here, 
requires  only  that  F  have  a  local  minimum  at  x°.  However, 
it  has  been  shown  ([23,  Theorem  8.13,  p.  368)  that  F(x) 
is  convex  on  its  domain — the  nonnegative  orthant— so  that 
every  local  minimum  is  in  fact  a  global  minimum,  and  the 
underlined  phrase  above  can  be  eliminated. 
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6.  Relation  of  Gibb s  Function  to  the  Hass  Action  Laws 

Let  x  be  a  c  ^position  vector  of  the  system  satisfying 
the  mass  balance  laws  (11.2.1)  and  the  nonnegatiyity  condi¬ 
tions  (11*2.2).  Such  a  composition  vector  we  shall  coll 
feasible.  We  wish  to  demonstrate  that  x  either  satisfies 
both  the  mass-action  laws  {(11.4.4)  or  (II.4.6))  and  the 
least-action  principle  for  the  free  energy,  or  that  x  sat¬ 
isfies  neither.  Actually^,  this  is  only  true  when  the  vec¬ 
tor  x  is  positive;  i.e.,  when  >  0,  j  *  1,  2,  ...,  n. 

If  for  some  j,  Xj  “0,  £hen  ‘the  mass-action  laws  no  longer 
make  sense  (except  in  trivial  cases) ,  so  that  no  comparison 
is  possible. 

This  question  is  dealt  with  at  length  in  [2 ] ;  we  will 
treat  here  only  the  simplest  case,  the  case  when  x  >  0, 

It  can  be  shown  that  £  is  minimized  .at  subject  to 
the  constraints  (11.2,1)  and  (11,2.2)  if  and  only  if  there 

exists  an  m-vceZor  tt°  such  that  x°  minimizes  the  Lngrangian 
function, 


(II.  6 d)  L(x,  tt°)  »  F(x)  -  tf°(Ax  —  b) 


subject  only  to  the  nonnegativity  conditions  (11.2.2/ 

But  x°  purely  cannot  minimize  L(x,  rr°)  subject  to 
x  >  0  unless  for  every  n-veotor  6  satisfying  x°  +  to  >  0 


for  t  >  0  sufficiently 
at  x°  in  the  direction 


small,  the  derivative  of  L(x,  t(°) 
Ci  (written  L*(x,  n°))  satisfies 


•»  JaHpr*-  __  ^(C>  ^  ^  v  *'  '  -  ~  ^  ^  '  jy  * 


(XX, 6.2)  Lg(x0,*-  -it0)  >  o. 

If  x°  >  0,  then  the  set  of  6  we  must  consider  includes 
every  n-yeetpr,  Farther,  it  can  be  shown  ([2],  Theorem 
8.11,  p.  368)  that  for  x°  >  0, 

(Ti.6.3)  L*<x°,  tt°)  -  |  Q^cj  +  log  *®  -  A*tt°) 

J  ^ 

where  Aj  denotes  the  column  of  the  matrix  A. 

.  Further,  from  (II.  6. 3)  we  see  that 

j  , 

(11.6,4')  ^(x0,  it0)  -  -Lj(x6,  it0) 

so  that  if  (II. 6. 2)  is  to  be  satisfied  both  for  0  and  —0, 
we  must  have  (if  x°  >  0): 

(11. 6. 5)  l’(x°,  n°)  =  0  V  e  e  En  . 

But  we  see  from  (II. 6. 3)  that  (II. 6. 5)  is  satisfied  for 
all  n-vcctors  0  if  and  only  if : 

(11. 6. 6)  c  +  log  Si0  «  ATrr°  . 

Compare  the  result  (II.6.6)  with  (II. 4. 5).  We  have  shown 
that  x°  >  0  satisfies  the  mass  action  laws  if  and  only  if 
x  minimizes  the  Gibb*s  function,  in  each  case  subject  to 
the  constraints  (II. 2. 3.)  and  (II. 2. 2). 
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7.  Nonide.il  Chemical  Systems 

In  nonideal  chemical  systems,  only  the  mass -action 
laws  change.  The  mass,  balance  laws  remain  linear.  It  is 
usual,  however,  to  cast  the  mass -action  laws  into  the  ideal 
form  by  substituting  functions  called  "activities"  for  con¬ 
centrations.  The  "activity"  of  a  species  j  is  thought  of 
as  its  "effective  concentration",  when  the  system  has 
composition  x,  and  is  a  function  of  the  composition.  Thus: 

(11. 7.1)  a^  *  a^Cx). 

As  will-  be  shown  in  Section  III. 7,  the  activity  is  a 
function  homogeneous  of  degree  zero  in  the  composition, 
so  that: 

(11. 7.2)  a..(x)  “  t  >  0. 

Often,  activity  coefficients  are  used.  The  activity 
coefficient  of  species  j  is  simply  the  ratio  between  its 
activity  and  its  concentration. 

a.  (x) 

(11. 7. 3)  \.(x)  =  - - 

J  X  • 

3 

It  can* be  seen  by  (II. 7.2)  and  (II. 1.2),  that  X.(x)  is 

3 

also  homogeneous  of  degree  zero. 

Using  activities,  the  nonideal  mass -action  laws  may 
be  stated  as  follows:  For  any  reaction  having  reaction 
vector  0,  there  is  an  equilibrium  constant  k(0),  such  th.  t, 
for  any  composition  vector  x  representing  an  equilibrium 
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state  of  the  system,  we  have; 

(11. 7.4)  n  (a.(x))6j  -  k(8). 

j-1  1 

The  same  development  as  used  in  Section  11.4  can  be 
used  to  show  that  there  must  exist  cons tan's  c?  such  that 

.  J 

if  x  is  an  equilibrium  compositionj  then  there  must  exist 
a  vector  n;  (n^,  . ..,  nn)  satisfying: 

(II. 7 *A)  c°  +  log  a(x)  »  A^tt. 

Alternately,  if: 

(11.7.5)  Cj(x)  »  c°  +  log  Xj(x), 
we  may  write  (II. 7.4)  as: 

(11. 7. 6)  c(x)  +  log  x  *  ATtt. 

It  can  be  shown  by  the  use  of  (II. 7. 2)  and  Euler's 
Theorem  on  homogeneous  functions  (see,  for  example,  [8] 
p.  234)  that  if  x  is  strictly  positive — i.e.,  x  >  0 — then 
x  can  satisfy  (II. 7. 4)  or  (II. 7. 6)  if  and  only  if  the 
function  F(x)  achieves  a  minimum  at  x  subject  to  the  mass 
balance  laws  (II. 2.1),  where: 

(11. 7. 7)  F(x)  »  S  x.(c?  +  loS  a. (x)) 

j=l  J  J  J 

n 

a  ^  xj(Cj(x)  +  log  X.). 

One  might  argue  that  some  x  satisfying  (11.7.4)  (and 
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(lit. 7,(5})  might  maximize  F(x)  rather  than  minimize  it. 
As  will  be  shown  later,  however,  F(x)  must  be  a  convex 

t 

function,  ruling  out  this  possibility. 


III.  THERMODYNAMICS 

The  Gibbs  function  Is  an  example  of  what  is  known  in 
thermodynamics  as  a  characteristic  function.  Characteristic 
functions  are  the  central  element  in  thermodynamics.  From 
the  characteristic  function  of  a  thermodynamic  system,  one 
may  compute  every  quantity  of  interest  to  the  thermodynam- 
icist.  In  this  section,  therefore,  we  will  explore  the 
concept  of  a  characteristic  function,  and  the  consequences 
of  the  laws  of  thermodynamics  on  such  functions. 

1.  Thermodynamic  Proper tie s 

Properties  of  a  thcrmodynaiv  .c  system  form  the  starting 
point  of  any  study  of  the  system.  They  are  defined  (more 
cr  less)  as  the  results  of  certain  measurements  done  on 
the  system,  and  certain  computations  done  with  these  mea¬ 
surements.  Thermodynamics  is  the  study  of  the  relation¬ 
ships  of  such  properties,  although,  strictly  speaking,  at 
least  one  property  of  interest  should  depend  on  whether 
the  system  is  hot  or  cold. 

It  is  not  our  purpose  to  define  these  properties  or 
to  describe  their  physical  interpretations.  We  arc  con¬ 
tent  to  list  a  few  of  the  more  common  ones.  These  include 
composition,  volume,  pressure,  energy,  temperature,  entropy, 
enthalpy,  chemical  potential. 
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2.  The  State  of  a  System 

Once  we  have  listed  all  the  properties  of  interest 
of  a  thermodynamic  system,  we  note  that  the  various  rela¬ 
tions  among  them,  as  described  by  thermodynamics,  allow 
some  to  be  computed  from  the  others.  Suppose  the  exper— 
imentor  chooses  from  this  set  of  properties  T  a  subset  S 
satisfying  the  two  conditions  below: 

(i)  Independence  —  it  should  not  be  possible 
to  compute  any  one  of  the  properties  in 
S  from  the  others  in  S. 

(ii)  Completeness  —  it  should  be  possible  to 
derive  every  property  in  T  from  a  know¬ 
ledge  of  only  thos  in  S. 

Then  S  is  one  possible  representation  of  the  state  of  Che 
system.  To  describe  the  state  of  the  syrtemn  one  must 
evaluate  each  of  the  properties  in  S. 

From  a  mathematical  point  of  view,  i  the.  properties  are 
simply  variables,  and  the  state  of  the  sysf,  ..m  is  described 
by  specifying  a  sufficient  number  (cc-mp J.c tones.;;  oi  ihu  - 
pendent  variables. 

Which  of  the  subsets  S  satisfying  (i)  and  (ii)  is  to 
be  chosen  is  a  matter  for  the  thermodynamicist  to 
decide,  on  tl  basis  of  convenience  and  good  taste.  Or¬ 
dinarily  he  would  choose  as  independent  those  variables 
most  easily  mea  arable  and  controllable,  so  long  ns  this 
choice  docs  not  violate  (i)  and  (ii). 
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For  example,  for  a  chemical  system,  most  of  whose 
phases  are  liquid  and  hence  incomprcss*.  blc,  one  would 
not  choose  volume  as  an  independent  variable.  One  would, 
rather,  choose  to  measure  and  control  the  pressure  of  each 
phase.  Similarly,  temperature  can  be  controlled  and  mea¬ 
sured  more  easily  than  internal  energy;  hence  one  chooses 
the  temperature  of  each  phase  as  another  independent  variable. 
The  description  is  completed  by  specifying  the  composition 
of  each  phase. 

In  more  complex  sy  terns,  a  state  description  might 
require  additional  properties.  For  a  more  complete  dis¬ 
cussion  of  these  matters,  sec  [ 1 ] ,  [6]. 


?? 


3.  Characteristi.c  Functions 


Suppose  we  have  listed  for  our  system  all  the  proper¬ 
ties  in  which  we  are  interested,  and  have  somehow  chosen 
from  among  them  a  subset  to  use  as  our  description  of  state. 
VJe  must  now  somehow  relate  the  state  variables  (i.e.,  the 
properties  used  in  the  description  of  the  state)  to  all 
the  others.  This  is  the  duty  of  the  characteristic  function, 
a  function  of  the  state  of  the  system  from  which  it  is 
possible,  by  various  manipulations,  to  compute  all  the 
properties  of  interest  which  are  not  themselves  state 
variables . 

In  16],  p.  24,  one  may  see  an  example  of  this.  Gug¬ 
genheim  shows  there  how  to  compute  from  the  Gibbs  function, 
its  tier  i\  a  elves  of  various  orders,  and  the  state  variables 
(in  this  case  temperature,  pressure,  and  composition  of 
each  i  base) ,  such  quantities  as  entropy,  enthalpy,  volume, 
energy,  chemical  potent!  Is,  and  others. 

Every  characteristic  function  of  a  given  system  is 
equivalent  to  every  other;  given  any  ch  racteriotic  func¬ 
tion,  explicitly  a  function  of  one  description  of  state 
Sp  we  arc  told  by  the  implicit  function  theorem  that  it 
can  also  be  considered  as  a  l unction  of  any  other  descrip¬ 
tion  of  state  $2 .  However,  for  particular  systems  and 
particular  problems,  one  characteristic  function  is  usually 
more  conv  nient  than  any  other.  In  our  casi ,  the  Gibbs  free 
energy  function  has  be*  a  found  most  appropriof  ■*. 
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Characteristic  functions  generally  have  the  inter¬ 
pretation  of  energy,  or  something  closely  related  such 
as  energy  divided  by  temperature. 
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4,  The  Global  Leas t-Ac l:lon  Principle 

One  of  the  uses  of  a  characteristic  function  is  to 
indicate  what  changes  in  a  thermodynamic  system  may  occur, 
and  what  changes  may  not.  This  is  because  if  the  thermo— 
dynamicist  chooses  wisely,  his  characteristic  function 
will  obey  a  global  least-action  principle. 

Given  a  then  odynamic  system,  one  may  observe  that 
its  state  changes  over  time.  It  may  be  in  state  s^  at 
time  tp  but  in  state  S2  t  at  time  t2«  A  function  F 
of  the  state  of  the  system  is  said  to  obey  the  global 
least-action  principle  if  for  any  two  such  states  and  S2, 

F(x2)  <  F(Sl). 

Thus  if  we  consider  as  the  system  evolves,  that  its  statu 
is  a  function  of  time,  then  F,  also  considered  a  function 
of  time,  must  be  monotone  decreasing.  (This  can  be  stated 
as  well  with  F  a  monotone  increasing  function  of  time. 

Then,  of  course,  — F  will  be  monotone  decreasing.) 

There  is  one  other  condition  that  F  must  satisfy  if 
it  is  to  obey  the  global  least— action  principle.  This  is 
the  local  least-action  principle  introduced  in  Section  IX, 
Thus  a  state  s°  may  be  a  stationary  state  if  and  only  if 
F  achieves  a  local  minimum  at  s°. 

In  this  section  v.c  will  show  that  every  characteristic, 
function,  i  f  chosen  wisely  for  the  syr.it..!  being  considered. 


obeys  the  global  least-action  principle.  It  is  only  neces¬ 
sary  to  assume  that  the  system  obeys  the  second  law  of 
thermodynamics.  We  shall  do  so  only  for  systems  with  one 
phase,  and  only  for  the  four  most  widely  known  and  used 
characteristic  functions.  In  fact,  it  would  be  sufficient 
to  demonstrate  this  for  a  single  characteristic  function, 
and  then  note  that  all  characteristic  functions  arc 
equivalent. 

We  will  require  the  following  thermodynamic  quantities: 

P  is  pressure 

T  is  temperature 

|i.  is  the  chemical  potential  of  species  j, 

J  a  function  of  the  state  of  the  system. 

V  is  volume 

q  is  heat  absorbed  by  the  system 

w  is  work  doi  •>.  on  the  system 
x.  is  the  amount  of  species  j, 

J 

The  characteristic  functions  we  will  consider  are: 

G  is  the  Gibb's  function 
II  is  enthalpy 

U  is  energy 

A  is  the  Kc1:i.holz  function. 


All  the  independent  infinitesimal  processes 
that  might  conceivably  take  place  may  be  divided 
into  three  types:  natural  processes,  unnatural 
processes,  and  reversible"  processes . 

Natural  processes  are  all  such  as  actually 
do  occur;  they  proceed  in  a  direction  towards 
equilibrium. 

An  unnatural  process  is  one  in  a  direction 
away  from  equilibrium ;  such  a  process  never 
occurs . 

As  a  limiting  case  between  natural  and 
unnatural  processes,  we  have  reversible  pro¬ 
cesses,  which  consist  of  the  passsagtTin 
eitTIef  direction  through  a  continuous  scries 
of  equilibrium  states.  Reversible  processes 
do  not  actually  occur,  but  in  whichever  direc¬ 
tion  we  contemplate  a  reversible  process  we 
can  by  a  small  change  in  the  conditions  produce 
a  natural  process  differing  as  little  as  we 
choose  from  the  reversible  process  contemplated. 


When  an  infinitesimal  process  occurs  in  a  thermodynamic 
system,  we  are  told  ([6],  p.  10)  that  the  change  dU  in  the 
energy  of  the  system  is  equal  to  the  sum  of  the.  ho. at  absorbed 
by  the  system  and  the  work  done  on  it.  That:  is, 

2' 


”v  - 


(111.4.1)  dU-q+w* 

If  heat  Is  lost,  or  if  work  is  done  by  the  system,  then 
q  or  w  appears  as  a  negative  number. 

The  Work  term  w  may  be  related  to  the  pressure  and 
volume  of  the  system  ([6],  p.  16) .  Thus: 

* 

(111. 4. 2)  w  «  -PdV. 


If  the  system  expands  by  an  amount  dV,  its  pressure  being 
P,  it  has  done  work  on  its  surroundings. 

The  Second  Law  of  thermodynamics  is  stated  in  terms 
of  entropy.  It  says  that  any  naturally  occuring  process 
will  increase  the  entropy  of  an  isolated  system.  An 
isolated  system  is  one  which  neither  mass  nor  energy  (in 
the  form  of  heat  or  work)  is  permitted  to  enter  or  leave 
the  system. 

When  the  system  is  not  isolated  the  change  in  entropy 
is  partitioned  into  an  external  and  an  internal  contribu¬ 
tion.  The  external  contribution  is  the  ratio  of  the  heat 
absorbed  to  the  termperature  at  which  the  heat  is  absorbed. 
Tine  internal  contribution  is,  according  to  the  second  law, 
nonnegative.  In  general,  then,  the  second  law  says  that 
for  any  natural  infinitesimal  process, 

(III. 4, 3)  dS  >  q/T. 
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Equality  holds  for  reversible  processes. 

The  entropy  function-  may  be  considered  as  a  charac¬ 
teristic  function,  in  fact  the  most  convenient  character¬ 
istic  function  to  u$e  for  systems  whose  energy,  volume  and 
material  content  are  constant.  Note  that  this  describes 
an  isolated  system,  for  which  the  second  law  takes  its 
simplest  form.  In  general  given  any  infinitesimal  process, 
the  change  in  entropy  may  be  computed  ([6],  p.  23)  by  the 
formula : 


(III. 4.4) 


T^dU  +  T^PdV  —  T”1  2  n,dx,. 

j  j  j 


Substituti  3  Equations  (III. 4.1)  and  (III. 4. 2)  into  (III. 4.4), 
we  see  that: 


(III. 4.5)  dS  =  S  -  l'"1  £  |i,dx,. 

I  j  J  J 

Temperature  in  thermodynamics  is  absolute  temperature. 
It  must  alv7ays  be  a  positive  number..  Thus  the  second  law 
simply  states  that  for  any  natural  process, 


(IIX.4.6)  S  | i.d:.  <  0. 

j  3  3  ~ 

Equality  holds  for  reversible.  processes. 
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As  previously  mentioned,  each  characteristic  function 
is  particularly  v;ell  suited  for  one  particular  description 
of  the  state  of  the  system.  If  the  state  of  the  system 
is  described  by  the  pressure,  the  temperature  and  the 
composition,  then  the  Gibb’s  function  is  best.  Then  for 
any  process, 

(111. 4.7)  dG  «  — SdT  +  VdP  +  2  p.dx.. 

j  J  J 

If  the  state  is  described  by  entropy,  pressure  and  com¬ 
position,  then  the  enth.'lpy  is  best: 

(111.4.8)  dH  «  TdS  +  VdP  +  S  [i-dx.. 

j  3  3 

If  the  state  is  described  by  entropy,  volume  and  composi¬ 
tion,  the  energy  is  best: 

(111. 4. 9)  dU  =  TdS  —  PdV  +  2  p.dx. 

j  3  3 

Finally,  if  the  state  is  described  by  temperature,  volume 
and  composition,  then  the  Hclmholz  function  is  best. 

(111. 4. 10)  dA  «  -SdT  -  PdV  +  2  p  .dx. . 

j  J  J 

Just  ns  the  second  law  stated  in  terms  of  entropy 
takes  a  particularly  simple  form  in  the  ease  of  an  isolated 
system,  so  if  it  is  stated  in  terms  of  different  characteristic 
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functions  will  it  take  a  simple  form  for  other,  nonisolated 
systems.  Thus  for  a  system  whose  temperature  and  pressure 
are  constant,  we  look  to  the  Gibb  *  s  function.  We  have 
dT  “  dP  -  0,  so  that  from  (III.4.7),  dG  becomes: 

(111. 4. 11)  dG»2|i.dx.. 

j  J  J 

If  entropy  and  pressure  are  constant,  we  have  dS  «  dp  «=  0, 
and  the  change  in  enthalpy  (III. 4. 8)  becomes: 

(111. 4. 12)  dH  »  51  HjdXj. 

3 

If  entropy  and  volume  are  constant,  then  dS  «  dV  =  0,  and 
the  change  in:  energy  (III. 4.9)  is  just: 

(111. 4. 13)  dU  «  2  u.dx.. 

j  3  3 

And  if  temperature  and  volume  are  constant,  then  dT  «  dV  =  0, 
and  the  Helmhclz  function  (III. 4. 11)  changes  by: 

(111. 4. 14)  dA  »  S  |i,dx,. 

j  3  3 

If  we  appeal  to  Eqiu  f:ion  (III. 4.6),  we  see  that  the 
second  law  may  be  restated  in  a  simple  form  for  any  system. 
Given  any  natural  process, 
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(a)  if  the  temperature  and  .pressure  are  constant 
throughout,  then  the  Gibb’s  function  is  monoton- 
ically  decreasing — i.e.,  obeys  the  least-action 
principle; 

(b)  if  the  entropy  and  pressure  are  constant  through¬ 
out,  then  enthalpy  obeys  the  least-action  principle; 

(c)  if  entropy  and  volume  are  constant  throughout, 
then  the  energy  obeys  the  least-action  principle; 

(d)  if  temperature  and  volume  are  constant  throughout, 
then  the  Holmholz  function  obeys  the  least-action 
principle. 
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5.  The  Form  of  Characteristic  Functions 

We  may  consider  a  thermodynamic  system,  to  be  completely 
defined  by  its  characteristic  function  F  (which  obeys  the 
least-action  principle),  and  a  set  of  constraints  on  the 
state  variables.  In  the  case  of  chemical  system,  the 
characteristic  function  is  the<  Gibb's  function,  and  the 
constraints  are  (1)  the  mass— balance  laws  and  the  non- 
negativity  conditions  on  the  composition  variables,  and 
(2)  the  constraints  that  pressure  and  tempera ture  be 
constant. 

If  this  is  all  we  know  about  our  system,  we  cannot 
be  content  with  Guggenheim's  classification  of  infinit¬ 
esimal  processes.  Maintaining  those  definitions  as  far 
as  possible,  we  say: 

Definition  III. 5.1:  An  (infinitesimal)  process  is  a  vector 
6  and  an  ini' ial  state  s°  satisfying  the  constraints  of 
the  system  such  that  s°  +  fcG  also  satisfies  the  constraints 
for  every  sufficiently  small  positive  real  t. 

In  our  case,  the  state  s  includes  composition  x  and 
temperatures  Ta  and  pressures  Va  of  each  phase.  Since  the 
Ta,  Pa  are  constrained  to  be  given  constants,  a  process 
0  must  have  its  changes  in  these  variables  equal  to  zero. 
Hence  only  composition  need  be  considered,  and  a  process 
reduces  to  a  reaction  vc  tor  such  that  x  +  to  >  0  for  t  >  0 
sufficiently  small.  Such  a  0  is  called  admissible  at  the 
composition  x. 
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Since  the  characteristic  function  satisfies  the 
least-action  principle,  we  know  that  any  process  that 
does  occur  leads  to  a  decrease  in  the  function;  Turning 
this  around,  we  define: 

Definition  III. 5. 2:  A  process  (s°,  0)  is  called  an  n-process 
(a  mnemonic  for  natural,  but  a  different  name  to  avoid 
confusion)  if: 

(III. 5.1)  litn  ns°)  <  o. 

t-0+  L 

It  is  called  a  u-process  (for  unnatural)  if: 

(III.  5. 2)  lim  Ehll-.tp  “ 
t-0+  C 

It  is  called  an  reprocess  (for  reversible)  if: 

(III. 5.3)  lim  -  F(s°)  =  o 

This  definition  limits  us  to  those  systems  whose  char¬ 
acteristic  functions  possess  all  the  limits  necessary  in 
the  definition;  but  this  includes  all  the  systems  that  we 
have  ever  seen. 

The  correspondence  between  Guggenheim's  processes 
and  ours  is  as  follows : 

(a)  A  process  is  unnatural  according  to  Guggenheim 
if  and  only  if  it  is  a  u-process. 

(b)  Every  natural  process  is  an  n-process. 

(c)  Every  reversible  process  is  an  r-proccss. 

The  converses  of  statements  (b)  and  (c)  are  not  true. 

Not  every  process  is  classified  by  Guggenheim. 
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Guggenheim  ([6],  p.  9)  states  the  first  law  as  follows: 

When  several  systems  interact  in  any  way  with 
one  another,  the  whole  set  of  systems  being  isolated 
from  the  rest  of  the  universe,  the  sum  of  the 
energies  of  the  several  systems  remains  constant. 

Here  the  word  "isolated”  means  that  neither  mass  nor  energy 
is  permitted  to  enter  or  leave'  the  set  of  systems. 

Suppose,  however,  that  we  permit  energy  to  enter  or 
leave;  the  set  of  systems,  but  close  the  systems  from  the 
entry  or  exit  of  mass.  Then  clearly  the  first  law  will 
still  hold  for  those  processes  such  that  no  energy  actually 
does  enter  or  leave.  We  have  carefully  defined  r— processes 
to  satisfy  exactly  this  condition. 

The  second  law,  of  course,  is  simply  the  statement 
that  the  characteristic  function  satisfies  the  least-action 
principle  for  the  system  under  consideration.  These  two 
laws  have  important  consequences  for  the  form  of  the  charac¬ 
teristic  function. 

Consider  a  chemical  system  consisting  of  a  single 
phase,  with  composition  x  and  characteristic  function  G(x) . 
We  also  call  G(x)  the  free  energy. 


If  wc  divide  the  phase  into  n  identical  parts,  then  the 
free  energy  of  each  part  is  the  same. 
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If  we  again  put  these  parts  together.,  we  find  the 
result  is  -indistinguishable  from  the  original  phase.  What 
we  have  done  and  undone  rrust  have  been  an  r-process,  so 
the  .first  law  applies.  That  is,  the  free  energy  of  the 
system  is  conserved,  and: 

G(-x)  -  ^G(x). 

Clearly,  G(^-  x)  is  the  free  energy  of  one  of  the  parts,  each 
of  which  has  composition  i  x. 

We  may  take  any  m  of  these  parts  (if  we  wish,  we  may 
let  m  >  n)  and  mix  them,  and  by  the  same  arguments  as  above, 

G(§  x)  -  §G(x). 


Characteristic  functions  are  always  considered  to  be 
continuous .  Thus  we  write  that  for  any  nonnegative  real 
number  a, 

(III. 5/4)  G(ax)  »  aG(x) . 

The  first  law  also  tells  us  something  about  multiphase 
systems.  Let  a  chemical  system  have  p  phases.  Let  the  kth 
phase  have  composition  vector  x^)  and  characteristic  (or 
free  energy)  function  G^(x^)  when  considered  as  a  separate 
system. 
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On  the  other  hand,  we  may  considi  >*  all  the  phases  to¬ 
gether,  as  though  they  are  a  single  system.  Iii  this  case 
the  characteristic  function  will  be  G(x^,  x^). 

But  such  a  'conceptual  change,  being  no  change  at  ail, 
is  an  r-process,  and  hence  conserves  free  energy.  Thus: 

(III. 5. 5)  G(x(1\  ...»  x<P>>  *  T.  G^(x(k^)  . 

k-1 


Property  1:  The  free-energy  function  (or  characteristic 
function)  of  any  chemical  system,  ideal  or  not,  can  bo 
partitioned  into  a  sum  of  functions,  one  for  each  phase  of 
the  system.  Each  of  these  subfunctions  is  a  function  only 
of  the  state  of  the  phase  for  which  it  is  written,  and  is 
homogeneous  of  df.„.*ee  one  in  the  composition  of  that  phase. 


Lemma  1 1 1. 5. 3:  The  ideal  Gibbs  free  energy  function, 
F(x)  =  >:  •  (c  ri-  log  S),  has  property  1. 


Proof;  In  this,  case,  if  F(x)  ~  E  p(^X(x(^)y,  we  see  that: 

k=l 

p(^)(x(^))  «  g  x.  (c.  +  log 

j  e<k>  J  1  J 


Multiplying  each  x,.  in  a  compartment  by  the  same  factor  a 
changes  not  at  all,  so  that  F^(x^)  is  clearly  homo¬ 
geneous  of  degree  one.  Q.E.D. 
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Next  consider  two  chemical  systems,  each  with  a  sing' e 
compartment.  'They  will  have  compositions  x  and  y  respectively. 


4 

We  can  arrange  their  free 

< 

energy  functions  in  a  suf- 

X 

, 

-4 

,  x  +  y 

ficieritly  general  way  to 

- 

:  t 

insure  that  the  first  has 

free  energy  G(x),  the 

second  G(y) .  Wo  can  do  this  by  allowing  in  G  for  an  energy 
contribution  from*  every  species  that  occurs  in  either 
compartment,  although  any  species  may  be  at  zero  level  in 
one  or  the  other  of  the  compartments. 

Experimental  evidence  indicates  that  if  x  is  not  a 
scalar  multiple  of  y,  then  mixing  the  two  compartments 
together  is  an  n-procoss.  In  this  event,  assuming  that 
x,  y  /  o, 

(III.  5. 6)  G(x  4*  y)  <  G(x)  +  G(y)  (if  x  r  ay)  . 

Definition  III.  5. 4;  Let  x  and  y  be  two  composition  vectors 
of  a  multicompar fr  onted  system.  Let  x^  be  the  composition 
of  the  compartment  when  the  whole  system  has  composition 
x,  and  similarly  for  y.  Then  x  and  y  are  said  to  be  quasi-. 
dependent  (written  x  ~  y)  if  for  each  compartment  <k>  there 
exist  numbers  a<k>,  p<k>,  not  both  zero,  such  that: 

a<lox^1C>  +  P<k>y(}°  =  0  for  every  compartment  <lo. 

This  notation  -was  first  introduced  in  [2],  p.  362. 
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Lemma  III. 5. 5:  If  x  ~  y,  then  for  every  species  j,  one  of 
the  following  holds: 

(i)  -  9j  (if  s<j>  >  o  and  y<j> >  0) 

(ii)  Either  x^  -  0  or  y^  *  0  or  both. 

Proof:  Clear  from  Definition  IV. 3.2.  Q.E.D. 

Using  the  notion  of  quasi-dependence,  and  inequality 
(III. 5. 6), coupled  with  the  Property  1.  in  a  simple  way,  we 
can  state: 

Property  2:  Let  x  and  y  be  two  compositions  of  a  multi- 
compartmented  chemical  system  with  free-energy  function  G. 
Then  G  is  linear  on  the  line  joining  x  and  y  if  x  ~  y, 
and  strictly  convex  otherwise. 

In  particular  (and  much  more  simply),  G  is  convex. 

Lemma  III. 5. 6:  F(x)  «  x  •  (c  +  log  &)  has  Property  2. 

Proof:  [2],  Theorem  8.13,  p.  368. 

It  should  be  pointed  out  here  that  the  author  realizes 
that  some  substances  are  only  partially  miscible.  An  example 
is  oil  and  water.  One  might  argue  that  such  a  system  is 
an  exception  to  Property  2,  since  a  phase  consisting  of 
oil  saturated  with  water  cannot  be  mixed  with  a  phase  con¬ 
sisting  of  wa tor  saturated  with  oil.  Indeed,  if  the  attempt 
is  made,  the  phases  will  spontaneously  separate.  Thus  the 
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free  energy  of  the  mixture  must  be  greater  than  the  sum 
of  the  free  energies  of  the  two  separate  solutions. 

This  proposed  counterexample  to  Property  2  neglects 
what  war  said  iii  Section  II.  1,  that,  molecules  of  the  same 
substance  but  in  different  phases  would  be  considered  as 
different  species.  This  is  not  simply  a  “cheap  way  out" 
of  the  dilemma i  A  molecule  in  a  solution  is  not  simply 
an  isolated  object.  It  has  interactions  with  other  mole¬ 
cules,  both  like  and  unlike,  in  the  solution.  It  is  these 
interactions  arid  the  consequent  configuration  of  the  poten¬ 
tial  field  surrounding  a  molecule,  as  well  as  the  atoms 
of  which  it  is  made,  that  determine  the  characteristics 
necessary  to  define  a  species.  All  those  characteristics, 
save  the  building  blocks  of  which  the  molecule  is  made,  are 
re-t resented  by  a  parameter  called  the  "chemical  potential," 
and  which  is  different  for  a  water  molecule  surrounded  by 
oil  molecules  than  for  a  water  molecule  surrounded  by  other 
water  molecules. 

More  detailed  information  on  this  poir.c  may  be  obtained 
from  any  book  on  statistical  mechanics. 

This  com] letes  the  exposition  of  the  desirable  properties 
of  characteristic  functions,  save  for  pointing  out  that  it  is 
always  assumed  that  whenever  x  >0,  they  are  (at  least)  twice 
continuously  differentiable,  and  that  they  are  continuous  on 
the  entire  sot  {x  |  x  >  0} . 


The  author  is  grateful  to  Dr.  N.  Z.  Shapiro  for  permission 


to  use  his  results,  the  derivations  of  Properties  1  and  2, 
in  this  papoi  . 
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6.  Boundary  Behavior  of  Free-Energy  Functions 

If  we  assume  that  x  >  0  and  compute  the  partial 
derivatives  of  the  ideal  Gibbs  free-energy  function,  we  find 
that. 


(III. 6.1)  jfc  -  c,  >  lP8  *4* 

j  J  J 

If  we  then  allow  x^  to  approach  zero,  while  retaining 
x^s,  >  0,  we  find: 

(HI. 6. 2)  Hm  ~  - 

xr°  ‘J 


It  appears  from  various  theoretical  considerations 
(e.g.,  [1],  p.  135)  that  (III. 6. 2)  is  generally  true  for 
arbitrary  free-energy  functions. 

We  have  explored  this  question  elsewhere  from  the  point  of 
view  of  statistical  mechanics,  and  it  appears  as  though  the 

only  assumption,  concerning  the  behavior  of  particles  in 
the  system  one  requires  in  order  to  prove  (III. 6. 2). is  that 
interactions  between  pairs  of  partic?  es  are  of  short  range 
only  (and  so  can  be  neglected  at  low  concentrations) . 

Let  x  be  a  composition  vector,  and  0  a  vector  such 
that  x  +  t0  >  0  for  sufficiently  small  positive  t.  One 
defines  for  such  pairs  x,  0,  and  for  each  free-energy  func¬ 
tion  F,  the  derivative  of  F  at  x  in  the  direction  0  to  be: 


(III. 6. 3)  F.(x)  «  lira,  - 

0  t-0 


F(x  -I-  te)  -  F(x) 


ia 


If  we  admit  extended'  real  numbers  then  this  limit  always 
exists  for  any  free-energy  function  corresponding  to  a 
real  physical  system.  (It  must  exist  for  almost  all  pairs 
(x,  0)  by  virtue  of  the  fact  that  F  is  convex. ) 

If  we  define: 

I  =  (<j>  i  x<d>  «  6} 

J  »  (j  I  x^  >  0  }: 

K  =  .{j  |  x^  >  0  but  Xj  =0}. 

then  since  free-energy  functions  are  always  highly  dif¬ 
ferentiable  and  because  of  the  homogeneity  property. 

(•III. 6. 4)  f'(x)  -  a  F<j>(e<j>)  +  a  o.  ,^<J><3>? 

9  <j>el  jcJ  J  Sxj 

if  e.  =  0  for  ovex-y  j.  e  IC>  and: 

3 

(III. 6. 5)  Fj(x).  =  -« 

if  for  some  j  c  K,  0.  >  0. 

It  can  be  shown  ([2],  p.  366),  that  the  ideal  free- 
energy  function 

F(x)  =  -xj(cj  +  log  *j) 
satisfies  (III. 6. 4-5). 

It  must  be  pointed  out  that  (III. 6. 4-5)  are  only 
assertions.  They  cannot  be  obtained  as  consequences  of 


previous  results. 
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7.  Implications  for  Chemical  Sys terns 

Consider  a  chemical  system  with  characteristic  func¬ 
tion  F(>:).  According  to  Sections  II. 5  and  II. 6,  this 
function  must  be  phase  separable,  and  the  function  written 
for  each  phase  must  be  homogeneous  of  degree  one,  and 
strictly  convex  except  on  rays  passing  through  the  origin. 


Let  F<k>(x<*<>)  be  the  function  associated  with  the 


th 


<k>™  phase.  By  Euler1  s  theorem  for  homogeneous  functions, 
(see  [8],  p.  234),  we  may  write,  for  x<k>  >  0  (the  interior 
of  the  domain  of  F^k>),, 


(III. 7.1)  F<k>(x<k>)  =  S  x. 


je<k>  I  " 


Summing  (111,7,1)  over  all  phases  we  find: 


(IT.I7.2) 


Ex 

j  ^ 


But  by  the  same  argument  as  produced  (III, 7.1)  we  know  that: 


(III. 7. 3) 


f(k)  =  S  x.  M 

J  CXj 


Clearly,  then,  In  addition,  since  F 
is  homogeneous  of  degree  one,  must  be  homogeneous  of 
degree  zero. 

For  each  j,  wc  arbitrarily  choose,  a  constant  c^, 
and  define: 
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(III. 7. 4) 


Substituting  (III. 7. 4)  into  (III.7.3)  shows  that  F(x)  may. 
be  written: 


(III. 7. 5)  F(x)  =  E.Xj(Cj  -I-  loga^x)). 

Compare  (III. 7. 5)  and  (II. 7. 7).  It  is  in  the  manner  above 
that  -the  activities  of  species  in  a  nonidcal  chemical  system 
may  :be  computed. 

We  'an  see,  though,  that  if  the  act  ivi tie's  are  computed 
.in  this  way,  they  must  satisfy  certain  conditions. 


(i)  From  Equation  (III. 7.4).  and  the  fact  that  F  is 

phase  separable,  we  find  that  a.(x)  is  homog  ncous 
of  degree  aci'o^  and  depends  only  on  the  variables 
x^  such  that  k  e  <j>. 

(ii)  Since  F(x)  is  convex,  we  must  have 


(III. 7. 6) 


>  o. 


If  the  phase  <j>  contains  more  than  one  species, 
then  the  infinitesimal  change  implies  by  partial 
differentiation  is  in  a  direction  along  which  V 
is  strictly  convex.  Thus  the  inequality  <"111.7.6) 
holds  strictly,  and: 


(III. 7. 7) 


(iii) 


(III. 7. 8) 


(III. 7. 9) 


B2F 

»■ 

dx/ 


1  Sa.(x) 

3j(xx  >  °- 


Since  a^(x)  >  0  (it  is  an  exponential),  aj(x) 
must  be  a  strictly  increasing  function  of  Xj, 
other  variables  being  held  constant. 

If  phase  <j>  contains  more  than  one  species, 
then  Xj  may  tend  towards  zero  while  x..^  remains 
bounded  away  from  sfero.  Thus  from  (III.  6. 2)  we 
haye:  (assum H '  >  ip 


■lira  EM 

x.-O  0K;j 

J  J 


lim  (c,  +  log  a ,  (x) ) 

Xj-'O  J  *J 


Clearly,  then. 


lini  a  .,(x)  «  0 . 
x.-*0  J 
J 

Thus  the  activity,  or  effective  concentration, 
of  a  species  must  be  zero  when  none  of.  that 

species  is  present— not  a  particularly  surpris¬ 
ing  result. 
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IV.  KH'IETICS 


Every  chemical  system,  and  likewise  every  thermo¬ 
dynamic  system,  is  a.  kinetic  system,  whose  evolution  can 
be  described  by  differential  equations  called  kinetic  laws. 
It  is  reasonable  to  ask  what  kinetic  systems,  correspond  to 
thcymodynamic  or  chemical  systems.  The  method  we  will  use 
to  examine  this  question  will  be  to  construct  from  the 
kinetic  system  a  function  satisfying  the  least— action 
principle.  Then  if  the  function  also  happens  l:o  have  the 
appropriate  pro  pc  las,  we  will  know  the  kinetic  system 
may  be  considered  a  thermodynamic  one  also., 
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X,  Kinetics, and  the -Least  Action  Principle 

The  chemical  system  considered  in  this  section  is  a 
special  case  of  what,  are  called  primitive  systems  in  [9], 

It  is  a  two-ph'is'e  chemical  system  satisfying  the  assumptions 

1.  There  are  no  charged  species. 

2.  There*  is  only  one  possible  reaction;  the 
migration  of  molecules  of  one  particular 
species  from  one  ph.  se  to  the  other. 

Molecules. of  all  other  species  are  re¬ 
stricted  to  remain  in  the  phase  they 
initially  occupy. 

3.  The  species  in  each  phase  are  miscible  in 
all  proportions. 

We  name  the  two  phases  Phase  I  and  Phase  II,  The 
composition  vectors  of  the  t\  o  phases  will  be  x  - 
(xp  x2,  ....  xn)  and  y  «  (yp  y2 ,  ...,  yn)  of  Phase  I 
and  of  Phase  II  respectively.  We  consider  that  a  molecule 
of  species  j  xn  Phase  I  (with  corresponding  composition 
variable  x^)  is  chemically  Indistinguishable  from  a 
molecule  of  species  j  in  Phase  II  (with  composition 
variable  y^) .  The  only  difference  is  that  they  occupy 
different  phases.  Note  that  if  we  wish  some  nonpermcable 
substance  (i.e.,  a  substance  not  allowed  to  pass  from  one 
phase  to  the  other)  to  appear  in  (say)  Phase  I  only,  we 


may  set  the  amount  of  the  same  substance  in  Phase  II  to 
zero..  Thus  we  are  not  kept  from  considering  mixtures  of 
very  different  compositions  indeed. 


~PhasG~’i  Phase  IT 


Far  notationai  convenience,  we 
assume  that  it  is.  -species  1  in 
each  phase  that  can  migrate  to 
the  other  phase. 


Consider  this  primitive  system  as  a  kinetic  system, 
that  is,  with  its  evolution  described  by  differential 
equations.  Further,  suppose  that  the  functions  expressing 
the  unidirectional  flows  of  species  1  from  Phase  I  to 

Phase  II,  and  from  Phase  II  to  Phase  I  are  knqv/n  to  be 

\  ' 

Jj  and  Jjj  respectively,  and  that: 


(TV.l.Ia)  Jj  »  JO-O 

(IV.  1,1b)  JIX«J(y). 

Thus,  the  unidirectional  rate  of  flow  of  species  1  from 
one  phase  to  the  other  is  assum'd  to  depend  only  on  the 
state  of  that  phase  from  which  the  flow  originates j  and 
the  function  describing  the  unidirectional,  flow  is  the 
same  regardless  of  the  direction  of  flow. 

Next,  v;e  sensibly  assume  that: 

J(x)  >  0 

if  8 


(IV. i. 2) 


<r- « 


for  all  compositions  x  >  0  (it  may  be  zero  if  x^  =  0) . 

Finally,  we  assume  that  J  is  continuously  differenti¬ 
able  in  each  variable,  whenever  that  variable  is  positive, 
and  to  avoid  a  trivial  problem,  that  the  gradient  of  J  is 
never  zero. 

Our  aim  in  this  section  and  tlie  next  is  to  construct 
from  this  primitive  kinetic  system  a  function  for  which 
the  global  least-action  principle  holds-. 

We  may  write  the  differential  equations  describing 
the  evolution  of  this  simple  system  as  follows: 

(IV.  1.4)  =  J(y)  ,  J(x) 

dy,  dx. 

(IV. 1. 5)  =  J(x)  -  J(y)  -  -  . 

We  wish  to  d  d  a  function  G  corresponding  to  the  function 
J,  such  that  G  has  the  following  properties: 

1.  G  should  be  phase  separable.  That  is,  we  assume 

rC» 

that  G  may  be  written  as: 

(IV. 1,6)  G(x,  y)  =  G(x)  +  G(y)  . 

We  expect  the  symmetry  of  G  as  a  consequence  of  the  symmetry 
of  the  kinetic  system. 
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2.  We  assume  G  to  be  continuous^  for  x  >  0  and  twice 
continuously  differentiable  for  x  >  0.  Thus  we  assume 
that  if: 


(IV. 1.7) 


?(x) 


_  aG(x) 


ox. 


then  ?■  and  its  gradient,  are  defined  for  x  >  0. 

3.  We  assume  G  satisfies  the  global  least— action 
principle  with  respect  to  J.  Thus: 


(rv  18)  0-  =  <  0 

uv.-uo;  ^  ay~  cfiT  ^  3yL  Ht^  5  u* 

We  further  assume  that  G  satisfies  the  global  least- 
action  principle  with  respect  to  J  nontrivially.  Thus 
==  0  if  and  only  if  : 


dxj  dy^ 


*  — -  =  0. 


We  obtain  an  immediate  consequence  of  these  three  assump¬ 
tions  by  substituting  (IV. 1.4-5)  and  (IV. 1.7)  into  (IV. 1.8) 
After  some  Simple  manipulations  this  yields: 


AO. 


(XV.  1.9)  -  ~f  =  t-5(x)  -  §(y)][J(x)  -  J(y) ]  >  0. 


By  as.sumpti.oiij  wo  have  equality  if  and  only  if: 


J(x)  -  J(y)  =  0. 
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Suppose  the  states  of  the  two  phases  are  very  nearly 
the  same;  then  we  let: 

x  *  y  +  te. 

Substitute  these  expressions  into  (IV. 1.9),  divide  by  t  , 
and  take  the  limit  as  t  -•  0.  We  find,  then,  that: 

(xv. l.io)  eT(ve)T(vj)e  >  o.  * 

We  have  taken  0  to  be  a  column  vector.;  0-  is  its  transpose. 
The  gradient  vectors  v$  and  VJ  are  considered  to  be  row 
vectors,  and  in  (IV. 1.10)  they  are  evaluated  at  the  point  y. ■ 
Equation  (IV. 1.10)  is  true  for  every  vector  0,  so 
long  as  all  the  partials  used  exist.  (Certainly  so  long 
as  y  >  0.)  Thus  the  matrix  (v$)^(vj)  must  be  positive 
semi-definite. 

Alternately,  we  may  vievj  (IV.  1.10)  as  requiring  that 
the  product  of  two  linear  functions  of  0  be  nonnegative. 
Clearly,  this  condition  is  satisfied  if  and  only  if  v$  is 
proportional  to  VJ  with  a  nonnegative  constant  of  propor¬ 
tionality.  It  is  obvious  that  if  »  kvj  for  some  k  >  0, 
then : 

(iv. l.ii)  e'r(v*)T('W)o  =  u||(vj)e||2  >  o. 
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*& 


Y  •  • 

V  r  - 

I '  . 

}  o 

S'  ,  .  * 

t 

i  t 

»  ,  i  - 

.;  If  on  the  other  hand",  (v$)  ^  k(V'J),  then  since  (VJ)  */•  Q  > 

t  J  -  '  ,  ; 

i  and  <V5)  ^  O'  (else  v?  -  k(vj)  with  lc  -  0),  the  two  sets  ■’ 

‘  { 0  j  (vj)  0  =  -1}  and  {o|(v|)o  •-  +i}'  describe  hyperplancs , 

-  *  *  ( 

which,  by  assumption  caiinot  be  parallel.  At.  any  value  of  ‘ 

I'  ;  8  where  they  meet,  (IV.  5. 10).  is  violated, 

}  t  3  j 

i  Finally,  if  vr.  -  k(vj)  for' lc  <  0,  then  (IV.  1.11)  is 

i  .  } 

!  '  still  true  with  the  final  inequality  reversed.  By  the  j  ■, 

assumption  that  (vj)  0,  there  will  exist  0  such  that  j 

t 

:  _  ( Vj)  0 '  r  0.  Hence  we  may  easily  violate  (IV,  1.10). 

■  We  may  therefore  write : 


» 

i 

i 


(IV. 1.12) 


vv(y)  •«  lc(y)  *  VJ(v) 


where  lc(y)  is  a  nonnegativc  (scalar)  function.  (See.  for 

example,  [8]  p.  322.)  'To  satisfy  our’  condition  that  G  be 

twice  continuously  differentiable;,  we  will  insist  that 

k(y)  Lc  continuous;  sad  to  avoid  trivial  G  and  satisfy 

do  d>:i  dy-\ 

the  condition  that  ^  <  0  unless  -v~  ~  =  0,  we  insist 

that  k(y)  be  strictly  positive. 

It  is  convenient,  as  we  shall  see,  to  choose  as  our 
integrating  function  k(y) ,  the  function  k(y)  - 
Substituting  this  choice  into  (IV. 1.12)  and  integrating, 
wc  find  that: 
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(IV. 1.13) 


»  $(y)  -'C  +  log  J(y) 


where  c  is  a  constant  of  integration. 

We  cannot  for  general  J  explicitly  integrate  (IV. 1.13) 
to  find  an  equation  for  G. 


X 


2 i  Generalized  Mass-?Action  Laws 

Let  us  consider  a  kinetic  system  with  a  single  reaction. 


(IV. 2.1) 


<2  r.N.  2  piN,. 

j-1  J  J‘  3=1  J  j 


The  forward  rate  of  this  reaction  shall  be  J7(x),  and  the 
backward  rate  ^(x).  Thus  if  species  j  is  a  reactant 
and  not  a  product  (r j  >  0,  p^  **  0) ,  then  its  rate  of 
disappearance  is : 


(IV.  2. 2) 


88  ^4‘T^(x), 

Dis 


J 


Similarly,  its  rate  of‘  appearance  is: 


(IV. 2, 3) 


dx 


%c~)  «  r  „  JT(x). 


Ap 


If  a  species  j  is  a  product  but  not  a  reactant- 
(r^  =0,  Pj  >  0),  then  the  similar  quantities  are: 


dx 


dx. 


(IV. 2.4)  ((j’F’’)  13  P ^  J  C 5  (tr;  )  83  P-?*J  (x)  • 

uu  Dis  J  An  3 


Ap 


For  spot-!.,  which  are  both  products  and  reactants 
(r^  >  0,  p^  >0),  the  appropriate  expressions  are  sums 
of  the  above. 

As  in  the  previous  section,  let  G(x)  be  a  least-actioi 
function  for  this  kinetic  system.  Then  G(x)  must  satisfy, 


5** 


r 


(IV*2*^  ||  ■=  Hr Hz’  ~  0 

with  equality  if  and,  only  if  55  ma^  substitute' 

from  Equations  (IV; 2. 2-4)  into  (IV, 2. 5)  to  find: 

(IV. 2. 6)  I  ||-[.(r,J-(X)  +  p.  J+(x))  -  (r.J+(x)  +  p.J“(x))]‘  <  0. 
j  — 1  j  **  •*  J 

Collecting  terms,  we  find  that: 

(IV. 2. 7)  [2  (r.  -  p,)§g-](.r(x)  -  J+(x))  <  0. 

3-1  J  J  xj 

Equality  wiil  hold  if  and  only  if  J~(x)  =  J**(x) .  , 

Looking  back,  one  may  see  that  Equation  (IV. 2. 7^  is 
quite  similar  to  Equation  (IV. 1.9).  We  define  a  reaction 
vector  0  as  before  by: 

(IV. 2. 8)  6j  -  (rj  -  Pj),  1  <  j  <  n. 

Then  the  development:  of  the  last  section  shows  that,  if 
we  find  G  satisfying: 

(IV. 2. 9)  (2  0.  |L)  -  log  iix2 

3=1  J  °x3  J-(x) 

then  G  will  be  a  least— action  function  for  the  kinetic 
system. 
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In  fact,  this  equation  is  exactly  the  equivalent  of 
one  of  the  forms  of  the  mass— action  laws  developed  in 
Section  II. 4.  We  have  seen  (Section  II,  and  [2],  p.  360) 
that  the  mass— action  laws  can  be  stated  as  follows: 

(IV. 2. 10)  tt  x!J  =  k(e) 

j=l  2 


where  as  before  0  is  a  reaction  vector.  The  mass— action 
laws  are  such  that  when  they,  the  mass— balance  laws,  and 
the  nom.egativity  conditions  on  x  are  all  simultaneously 
satisfied,  then  the  system  is  in  equilibrium. 

These  laws  are  easily  derived  from  ideal  chemical 
kinetics  (and  from  many  other  forms  of  chemical  kinetics) 
as  follows. 

Fcr  the  reaction  (IV. 2.1),  the  simplest  form  of 
chemical  kinetics  (10,  p.  10)  assumes  that  the  unidirec¬ 
tional  rate  of  disappearance  of  reactan  •  j  can  be  written: 


(IV. 2. 11) 


and  that  its  unidirectional  rate  of  appearance  can  be 
written  as: 


(IV. 2. 12) 


(a r suming 


that  p. 
1  J 


0). 


5b 


The  numbers  P^,  P_-j_  are  called  respectively  the  forward 
and  backward  rate  constants  for  the  reaction.  Similar 
expressions  hold  for  those  x^  which  are  products. 

Taking  the  ratio  of  (IV.2.11)  to  (IV.2.12)  and  noting 
that  at  equilibrium  species  j  should  appear  exactly  as 
fast  as  it  disappears,  we  see  that  at  equilibrium. 


(IV. 2. 13) 


n 

n 

j=i 


xj 


(rrpj) 


-l 


If  we  let  8,  »  r.  -  p.  for  each  j  «  1,  2,  n  (as  in 

J  J  J  p 

Section  II)  and  k(e)  =  we  see  that  (IV. 2. 13)  is 

H1 

identical  to  (IV. 2. 10). 

Note  that  we  arrive  at  the  same  expression  even  if 
we  multiply  both  P_x  and  Px  by  an  arbitrary  positive 
function. 

On  the  other  hand,  if  the  system  is  not  at  equilibrium, 
it  is  still  true  that: 


n  n  q  P  -j 

(XV. 2. 14)  E  9,  =  log  (  n  X?)  pi- 

j=i  J  oxj  Lj=1  J  _1. 

where  F  is  the  original  Gibbs  function.  It  is  clear  that 
Equation  (IV. 2. 14)  is  the  same  as  Equation  (IV. 2. 9)  in  the 
case  of  ideal  chemical  kinetics,  and  that  the  analog  of 
the  mass— action  lav;  in  the  case  of  a  more  general  kinetic 
system  is  that  if  x  is  any  equilibrium  composition,  then: 


(IV. 2. 15) 


j*(x)  „  , 

jrw 
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3.  Kinetic  Systems  with  Several  Reactions 

In  general,  a  kinetic  system  will  have  more  than  one 
reaction.  Let  us  say  that  it  has  a  total  of  K  reactions, 
numbered  1  <  k  <  K,  each  one  of  the  form: 

(IV. 3.1)  S  rP^N.  ^  S  pP^N.. 

j=l  J  J  j-1  J  3 


Typically,  an  experimentor  will  not  measure  the  forward 
and  backward  rates  of  each  reaction.  Instead,  using  tracers, 
or  some  other  marking  method,  he  will  measure  the  unidirec¬ 
tional  rates  of  appearance  or  disappearance  of  various  species. 
From  these  quantities  we  must  compute  the  act-.;.;l  reaction 
rates. 

Clearly,  the  equations  relating  the  rates  of  appearance 
and  disappearance  of  the  species  and  the  reaction  rates  are 
sums  of  terms  of  the  sort  (IV. 2. 2-4).  Thus: 


(IV. 3. 2) 


(^AP  ■  }  '(k> 


K 


r^'JT(x)  +  E  p«j+ 
1  3  K  k-1  3  k 


(x) 


dx.  K 

%t^  **  s 

‘  at:  Dis  k 


.(k)  T+ 


K 


i  rV^jJ(x)  +  E  P<k>.T(x) 
1  J  k**!  J 


e 


We  assume  that  the  reaction  vectors  0  ^  (defined  by 

*  rj^  “  P t  1  <  j  <  n  are  linearly  independent. 
In  this  case,  whenever  Equations  (IV, 3.2)  possess  a  solu¬ 
tion,  it  will  be  unique.  Thus  if  we  let: 
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dx. 

(IV. 3.3) 


dx. 

(he1) 


Ap 


1  <  j  <  n 


then  we  may  write  the  difference  of  the  two  Equations 
(IV. 3. 2)  as: 

(IV. 3.4)  ^  -  s  ^k)fk(x)  1  <  J  <  “ 

KBJL  • 

where  ffc(x)  =  J£(x)  -  j£(x) .  If  (17.3.4)  has  a  solution 

(i.e.,  if  ^  is  a  linear  combination  of  the  reaction  vectors 

e(k)),  then  the  solution  is  unique.  A  properly  formulated 

kinetic  problem  will  surely  possess  a  solution.. 

It  is  a  simple  matter  to  substitute  the  solution  f(x) 

to  (IV. 3. 4)  into  either  of  Equations  (IV. 3.2)'.  *We  choos-. 

to  use  the  second  of  thc::e  equations.  Then: 

•  • 


(IV. 3. 5)  (£i)Dis  -  jL  r<k>£k(x)  -  ?f)^(x). 

Again,  if  Equations  (IV. 3. 2)  arise  from  an  actual  kinetic 
system,  (IV. 3. 5)  must  possess  a  solution  J£(x),  1 <  k  <  K. 

Once  we  have  arrived  at  expressions  for  the  unidirec¬ 
tional  reaction  rates  J+  and  J'”,  we  may  use  the  results 
of  the  last  section.  Thus  if  G  is  to  be  the  least— action 
function  for  this  kinetic  system,  we  require  that  G  simul¬ 
taneously  satisfy  the  K  independent  equations. 


(IV. 3. 6) 


n 


E 

j-1 


n(k)  H- 

j  ox . 

•  3 


log 


Jk(x) 


1  <  k  <  K. 
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4.  Kinetics  and  Thermodynamics 

Given  a  kinetic  system  such  as  that  described  in  the 
last  section,  we  wish  to  decide  if  it  may  be  considered  a 
thermodynamic  system.  To  be  so  considered,  the  function  G 
satisfying  conditions  (IV. 3. 6)  must  have  the  usual  thermo¬ 
dynamic  properties:  it  must  be  phase  separable  and  homo¬ 
geneous  of  degree  one;  the  function  G*'^  for  a  particular 
phase  <lo  must  be  strictly  convex  except  on  lines  through 

the  origin;  and  -  must  behave  in  the  appropriate  way 

3 

as  x.  approaches  zero.  These  conditions  on  G  will  imply 
conditions  tn  the  kinetic  laws  J  and 


Positivity:  For  Equation  (IV. 3. 6)  to  be  defined,  the 

ratio 

jjifo) 

must  be  positive.  This  will  be  assured  if  the  kinetic 
system  is  a  real  system,  for  in  that  case  wc  will  have  the 
following:  Let  the  raction  associated  with  reaction  lc  be: 


n  Pc') 

?!  E  p!j  N, 
j=l  3  3 


with  forward  reaction  rate  J^.(x),  and  backward  rate  J^(x)  . 
Then: 


(IV. 4,1)  4<x)  >  0  and  Jj(x)  >  0  vx, 

and 

(IV. 4. 2a)  j£(x)  «  0  if  and  only  if  xj  “  0 

for  some  j  such  that  rj^  >  0, 

(IV. 4. 2b)  J£(x)  »  0  if  and  only  if  xj  “  0 

for  some  j  such  that  >  0. 

Continuity  an'1  ^5  f ferret lability?  We  will  insist 
that  G  be  tv, 'lee  continuously  differentiable  for  x  >  0, 

immediately  from  (IV. 3. 6) 


for  x  >  0.  It  is  not 

must  be  continuous  for  x  >  0, 


and  continuous  for  x  >  0.  One  sees 
that  if  this  is  co.  then 

i«o  c-i— ) 

J^(k) 

must  be  continuously  differentiable 
true  that 


log  (-£ 
% 


6l 


Homogeneity  and  Phase  Scpavflbl3.lt/:  If  G  is  to  bo 

<)G 

homogeneous  of  degrees  one  and  phase  separable,  then 
must  be  homogeneous  of  degree  zero  and  phase  separable. 
From  (IV. 3.0)  we  sec  that  fur  the  kinetic  system  to  be 
considered  a  thermodynamic  system,  the  kinetic  lows  j£ 
and  must  satisfy  the  following  conditions:  Let  x  and 
y  satisfy: 


(IV. 4. 3) 


V<j>  “  V<j>'  1  *  J  s n- 


One  may  easily  verify  that  G  should  be  linear  on  the  line 

joining  x  and  y  Ilcnce  along  that  line,  must  be 

oxj 

constant.  Thus  (IV. 3. 6)  implies: 


(IV. 4. 4) 


,4<*K  A^k 

log  (•—  log  (——•—■)  . 

JfeOO  Jfe(y) 


Notice  that  neither  J*(<)  nor  J^(;:)  need  be  homo¬ 
geneous  of  any  degree.  The  restriction  is  on  their  ratio 
only. 

Convexity;  Lot  8  be  any  reaction  vector.  Clearly 

we  may  express  <)  as  a  linear  combination  of  the  vectors 
K 

(0V  '}  .  Wc  suppose  that: 

k«l 


(IV.4.3) 


n  rj  l  v.  o 
lo-l  lc 
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If  G  is  to  be  convex,  wo  know  that: 


2 

(IV. 4. 6)  3-2  G(x  +  80)  Is  -  0T(v2G(x))e  >  0 

ds 

where  v2G  is  the  matrix  of  second  partial  derivatives  of  G. 
Substituting  (IV. 4. 5)  into  (IV.4.6)  we  find  that: 


(IV. 4. 7) 


K 

E  V. 

k-1  K 


n 

E  0 
i-1 


(k)JL 

i 


n 

E 

>1 


e 


(4)  3G 
J  ex. 


The  quantity  in  square  brackets  is  recognizable  from 
(IV. 3. 6)  as 


log  ( 


Therefore,  if  we  define: 


(IV. 4. 8)  |ikt(x) 


E 

i-1 


.<10/1  i 

1  \jfoO  °’xi  ~J ZM  °\  I 


and  if  the  matrix  U(x)  is  taken  to  be  8  K  x  K  matrix  whose 

i.L 

element  is  \x^,  then  (IV.4.6)  simply  states  that 
U(x)  is  a  poritivo  sernidofinitc  matrix. 

Wc  also  knov;  that  G  is  to  be  strictly  co- vex  wherever 
phase-separability  and  homogeneity  do  not  imply  that  it 
is  linear.  Thus  for  any  vector  v^  (v^,  V£,  ...,  vk). 


(IV. 4. 9)  vT(U(x))v  >  0 
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whenever  0  *  E  v,.0^  fails  to  satisfy: 

k*i  K  - 

(IV. 4. 10)  0<j>xj  "  0jx<j>  1  <  J  < 

But  what  if  6  satisfies  (XV. 4. 10)?  In  that  case 


we  notice  that: 


n 


(XV.4.11)  ^  VW 


E  vk0 


00 


[7  log  ( 


jt(x) 


•  tfs  ^log 


■fo*  +  s0\  . 

J7(x  *1*  so)/ 


)3 


s«*0 


due  to  condition  (IV. 4.4).  lienee  for  those  v  we  must  have 
v  ll(x) v  >o  0j  as  required. 

It  can  easily  be  verified  that  in  the  case  of  the 
primitive  kinetic  system  of  Section  IV, 1,  Equation  (XV.4. 6) 
reduces  to  the  condition  that  the  rote  of  flow  of  a  species 
out  of  its  phase  must  he  a  monotone  increasing  function 
of  the  amount  of  that  species  present,  other  things  being 


equal. 

Notice  also  that  U(::)  must  be  a  symmetric  matrix. 

Thus: 


6U 


(XV.  4. 12)  \ikl 


£  o>K/  <-* 


-  E  6 


,0 0 

i-1  1  Jj(x)  ~Ai 

n  (k)  '  n 


aji 


ox. 


ajr 

~-i) 

ox. 

X 


9  (•  r  „(*)  ac 

i-i  1  ~*7  j-i  0J 


OX 


j 


I  6<*> 


j-1 

n 

£  0 

j-1 


’j 


JL  (  E  fi(k) 

0Jij  i»l  1  oxi 


'  t-tv.a  ox7  ™  oxT' 


J£(*) 


* 


-  M 


4k  ‘ 


Boundary  Behavior:  For  some  fixed  k,  let  Gjk^  >  0, 
(We  could  as  easily  use  some  ejk^  <  0.)  Let  x  be  such 
that  each  x^  is  greater  than  zero,  i  v4  j  and  tends 
toward  zero.  Then  v?c  must  have: 


(XV. 4. 13) 


n 

lim  F.  o> 
x.-O  i«l  L 

xi>0'J 


00  SG 
"6y~, 


Substituting  from  (XV. 3. 6)  v:e  find  that  the  kinetic  aysti 
cannot  be  a  thermadynatrie  system  unless: 

j/,(>0  for  ('S1^  >  0 

(XV.4.14)  lii:i  (log  J* — -)  »  \  4.  , 

"*<-0  ,  JjTC-O  for  rV:'  <0. 

-  >0  i/  j  K  J 


Consider  the  re  a ’t  I  on; 
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S  r/^N,  S  p/k^N.. 


1=1  J 


j=l 


Clearly,  for  any  species  j  such  that  r^')  >  o 

p  (h)  „  or  £or  mhich  r.^c^  »  0  and  p.^c^  >  0,  Equations 
J  J  J 

(IV.4.2a-b)  imply  (IV. 4. 14).  Only  in  the  case  that  both 

r.^  >  0  and  p.^k)  >0,  and  r/1^  £  p , ^  does  (IV. 4. 14) 

J  3  J  J 

constrain  the  kinetic  system* 

Cycling:  Suppose  there  is  an  x°  such,  that  if  the 
kinetic  system  is  initially  at  >:°  then  it  will  eventually 
return  to  >:°,  having  passed  through  states  y  !}:  x° .  This 
phenomenon  is  called  cycling.  A  kinetic  system  which  can 
cy<'lc  cannot  be  a  thermodynamic  system,  for  the  thermo¬ 
dynamic  function  C  would  have  to  decrease  strictly  through¬ 
out  the  cycle,  and  the  value  of  G  would  not  be  wall  defined 
at  any  point  y  in  the  cycle. 

Even  when  the  kinetic  laws  J1  and  S~  satisfy  all  the  . 
above  criteria — positivity,  differentiability,  ctco — it 
may  be  difficult  to  construct  a  function  G  for  the  thermo¬ 
dynamic  system.  Indeed,  we  have  not  established  that  for 
all  such  J  a  function  G  must  exist. 

V’hen  a  kinetic  system  does  admit  a  thermodynamic 
description,  however,  there  is  some  assurance  that  its 
least-action  function  may  be  interpreted  as  a  measure  of 


and 


s 

4 
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free,  or  usable,  energy.  One  can  imagine,  for  example,  a 
machine  whose  motive  power  is  obtained  from  a  reaction 

proceeding  in  one  direction  faster  than  the  other,  and 

i 

xdiich  yields  useful  work  in  the  process.  Thus  there  is 
some  physical,  as  well  as  a  mathematical  basis  for  calling 
a  kinetic  system  that  admits  this  type  of  description  a 
thermodynamic  system. 
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